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ABSTRACT 



Context. TMR-1 (IRAS 04361+2547) is a class I proto-stellar source located in the nearby Taurus star- forming region. Its circumstel- 
lar environment is characterized by extended dust emission with complex structures and conspicuous filaments. A faint companion, 
called TMR-1 C, located near the proto-star had been detected in previous studies, but its nature as a very young substellar object 
remained inconclusive. 

Aims. We aim at improving the constraints on the nature of the faint object TMR-1C, and to investigate the process of very low-mass 
star formation in the TMR-1 system. 

Methods. Using very sensitive infrared imaging observations of the TMR-1 system as well as near-infrared spectroscopy, we construct 
the spectral energy distribution of TMR-1C over a much larger wavelength range as had been possible in previous work. We then com- 
pare the spectral energy distribution with models of extincted background stars, young sub-stellar objects, and very low-mass stars 
with circumstellar disk and envelope emission. We also search for additional low-luminosity sources in the immediate environment 
of the TMR-1 proto-stellar source. Furthermore, we study the surrounding near-infrared dust morphology, and analyse the emission 
line spectrum of a filamentary structure in the physical context of a bow-shock model. 

Results. We find that the observed spectral energy distribution of TMR-1 C is inconsistent with an extincted background star, nor can 
be fitted with available models for a young extremely low-mass (< 12Mj) object. Our near-infrared spectrum indicates an effective 
temperature >3000K. Based on a good match of TMR-lC's spectral energy distribution with radiation transfer models of young stellar 
objects with circumstellar disks, we propose that TMR-1C is most likely a very low-mass star with Ms 0.1 - 0.2M o surrounded by 
a circumstellar disk with high inclination, i > 80°. Interestingly, we detect an additional very faint source, which we call TMR-1D, 
and that shows a quite striking symmetry in position with TMR-1C. TMR-1 C and TMR-1D may have been formed from a common 
triggered star-formation event, caused by a powerful outflow or by the collision of primordial proto-stellar disks. The impact of an 
outflow is traced by molecular hydrogen emission that we detect from a distinct filament pointing towards TMR-1C. A comparison 
with C-type bow shock models confirms that the emission is caused by shock excitation. 

Key words. Stars: low-mass, brown dwarfs - Stars: pre-main sequence - Shock waves - Stars: individual: IRAS 04361+2547 
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<0 1. Introduction on the basis of the presence of a striking filament structure 

^2 . that arises from TMR-1 and points directly towards TMR-1 C. 

^H . TMR-1 (IRAS 04361+2547) is a class I young stellar object lo- T98 further sugges ted that TMR1-C was catapulted to its cur- 

L; ; cated in the Taurus molecular cloud, that was actually resolved rent location due to dynamical interactions with the proto-binary 

. ■ r-j into a binary source with a measured components' separation TMR1-AB and that the arc-shaped filament could trace the ejec- 

^ . of 0.31"(Terebey et al. 1998, hereafter T98). We will therefore tion path of TMR . 1C through the gaseous Mailing circumstel- 

H ; refer to it as TMR-1AB in the following. The total bolometac lar envelope of TMR-1. The very low luminosity suggested for 

luminosity of TMR-1AB had been estimated to ~ 2.8L Q , indica- TMR _i C indicated that it should be a substellar object, maybe 

live for a low-mass protostellar system (Kenyon et al. 1993). The even a p lane tary mass object. However, the physical association 

circumstellar environment of TMR-1 AB is characterized by ex- of TMR . 1C with T MR-1, and hence its nature as a substellar 

tended emission from a dusty proto-stellar envelope and from ob j ecti was strong i y debated during the years after its discov- 

patches of molecular cloud material left over from the proto- ery In an attemp t to clarify the evolutionary status of TMR-1 C, 

stellar collapse. In their sensitive HST/NICMOS observations, Tere bey et al. (2000) carried out near-infrared spectroscopy us- 

T98 detected a faint compact object, named TMR-1C, at a pro- ing the Keck telescope. The result of these observations showed 

jected distance of ~ 10" from TMR-1AB, which corresponds to that the spe cti-um of TMR-1C, at the signal-to-noise level that 

-1400AU at the distance of the Taurus molecular cloud. The could be reac hed, is consistent with an extincted background 

physical association of TMR- 1 C with TMR- 1 AB was suggested dwarf star spe ctrum, but still room was left for an interpretation 

within the extremely low-mass object ejection hypothesis. 

* Based on observations made with ESO telescopes at the La Silla 
Paranal Observatory under programme ID 265.C-5747(A), and based In this paper we use ESO Very Large Telescope (VLT) data 

on data obtained from the ESO Science Archive Facility obtained with the Infrared Spectrometer and Array Camera 
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(ISAAC), as well as Spitzer/IRAC observations in order to re- 
visit the TMR-1 system. After the description of the observa- 
tions and the data reduction (Sec. 2), we discuss in Sec. 3.1 the 
morphology of the circumstellar dusty environment of TMR- 1 
and report on the detection of new features and objects identi- 
fied in our sensitive ISAAC images. We then use our K-band 
low-resolution spectroscopy together with the spectral energy 
distribution of TMR-1C, which we construct from the photom- 
etry presented in this paper and collected from the literature, in 
order to analyze the nature of TMR-1C in Sec. 3.2 and Sec. 3.3. 
Our K-band spectroscopy was also performed on a significant 
part of the filament structure "connecting" TMR-1C with TMR- 
1 AB. Since the spectral resolution is almost 4 times higher than 
in previous spectroscopic observations we are able to resolve nu- 
merous emission-line features arising from the filament. In Sec. 
3.4 we analyze the filament spectrum in detail, which leads us to 
conclude that a significant part of its emission is characterized 
by shock-excited emission. In Sec. 4 we present scenarios that 
could possibly explain the symmetry of the circumstellar nebu- 
losities, the filamentary structures, and the existence of a pair of 
very low-mass objects, as a physical entity. Finally, we summa- 
rize our results and conclusions in Sec. 5. 

2. Imaging and Spectroscopy Data 

2.1. Near-infrared photometry from VLT/ISAAC HK S 
observations 

A first set of images of TMR-1 was obtained at Ks-band (A c - 
2.16/an, AA = 0.27/zm) with ISAAC at the VLT-ANTU(UTl) 
telescope as part of commissioning of the instrument. These data 
were taken during the night December 04-05, 1998, immedi- 
ately before the spectroscopic observations, providing also a flux 
calibration for these (Sec. 12.3b . The data are accessible via the 
ESO Science Archive Facility^. The scale of the IS AAC/Hawaii 
Rockwell detector was 0.147"/pix, yielding a field-of-view of 
2.5'x2.5'. The images typically consisted of 50 frames averaged 
with 1.77 seconds individual exposure time. Those having the 
best image quality, of the order of ~ 0.55" FWHM, were se- 
lected for further data reduction. The sky contribution was deter- 
mined from median combining the individual images, and then 
subtracted from each image. Finally, the sky subtracted images 
were badpixel corrected, aligned and combined, and trimmed to 
a common overlap region. The total integration time in the final 
mosaicked image was 177 sec. 

A second set of ISAAC near-infrared data was obtained 
at H-band (A c = 1.65/im, AA - 0.30yum), and again at Ky- 
band, on two nights during October 2000 (TableQ]). The images 
were taken by using a random jitter pattern centered on TMR- 
1 . In order to avoid strong saturation from the brightest source, 
TMR-1 AB, short detector integration times of DIT=2sec and 
DIT=3sec were used for K$ and H exposures respectively. The 
total exposure times per image were 2x20(DITxNDIT)sec for 
Kg images and 3x20(DITxNDIT)sec for H images. The data 
were then processed with the Eclipse^ data reduction package 
(Devillard 1997), using the jitter algorithm. The images were 
dark-subtracted and flat-fielded with a skyflat obtained from ex- 
posures taken during twilight, then sky-subtracted, and finally 
recentered and stacked to single deep H- and Ks-band images. 
The total integration time in the stacked H-image is 12 minutes, 
while it is 20minutes for the Ks-image. TMR-1 and its imme- 
diate surrounding as seen in these sensitive images is shown 

1 http://archive.eso.org/eso/eso_archive_main.html 

2 http://www.eso.org/projects/aot/eclipse/ 



in Figures QJ, b, and c. The measured image qualities are very 
good, with FWHM- 0.5" (Ks) and FWHM- 0.7" (H), but still 
not good enough to resolve the binary TMR-1AB. 

Accurate photometry for all sources with a minimum S/N 
detection of 8 was then performed on the H and Ks mosaic im- 
ages by using PSF-fitting with IRAF/daophot. In particular for 
TMR-1C, which is close to the nebulous filament and the gen- 
eral nebulosity surrounding TMR1, PSF photometry provides 
a much more accurate photometric measurement than aperture 
photometry. In all of our final Ks and H mosaics TMR-1C has 
been well detected at a signal-to-noise ratio of >15. For the ob- 
servations carried out in 2000, the absolute photometric cali- 
bration has been derived from the UKIRT faint standard FS12 
that was observed immediately after our observations of TMR- 
1. The resulting photometry for TMR-1C is Ks=17.53mag and 
H=19.21 mag with an uncertainty of ~ 5%. We checked the pho- 
tometric calibration against the photometry of a 2MASS star that 
is present in the ISAAC field (2MASS 04391 199+2553490), and 
find agreement within 0.06 mag. Our Ks photometry of TMR-1C 
is also consistent with values reported by other authors at sim- 
ilar wavelength (Itoh et al. 1999), and also in agreement with 
the Ks-photometry determined from the 1998 imaging data (see 
Table [TJ. The photometric accuracy for the 1998 data is worse 
due to the lack of standard star observations during the observing 
night, so that absolute photometry was derived assuming the Ks- 
zeropoint determined for the ISAAC commissioning I period. 

Our deep imaging photometry revealed numerous faint ob- 
jects, detected all over the field-of-view (Figure [T^). There are 
more than 20 objects with a Ks-magnitude similar or fainter than 
TMR-1C, indicating that we are beginning to penetrate through 
the large amount of dust of the molecular cloud core. In this re- 
spect, the faintness and appearance of TMR-1C seems not to be 
any special, except for its location at the end of the filament. 
Very interestingly, we also detected a previously unknown, faint 
object which looks like a TMR-1C counterpart: it is located 
close to the end of a broad filamentary structure that extends 
from the large nebulosity seen to the North- West of TMR-1AB 
(FigureQJ, b and c). The object is clearly a > 10cr detection in 
the deep H and Ks-band images and fainter than TMR-1C by 
~ l.lmag at Ks and by ~ l.Omag at H. We will call this object 
TMR-1D hereafter. Although at this point we don't have a proof 
that there is any physical relation between this object and TMR- 1 
or to TMR-1C, we find it quite striking to detect such a symme- 
try. In the course of this paper, we will discuss the presence of 
TMR-1D in the context of different scenarios. 



2.2. Mid-infrared photometry from VLT/ISAAC L '-band and 
Spitzer/IRAC observations 

ISAAC mid-infrared imaging at L'-band (A c - 3.78/im, AA = 
0.58/im), had been carried out with the long-wavelength imag- 
ing camera of ISAAC, equipped with an Aladdin 1024x1 024 pix 
array, and providing a pixel scale of 0.071 "/pix in LWI4 mode. 
Data were taken in uncorrelated/high bias read-out mode using 
standard chopping and nodding techniques, with a chop-throw 
position angle of 122°, an amplitude of 30", and a chopping fre- 
quency of 0. 1 Hz. The reduction process included bad pixel fil- 
tering and flat-fielding. Then, positive and negative beams from 
the nod A and B positions were combined. The total integra- 
tion time yielded for the final image, shown in Figure Id, is 
~ 3000sec. TMR-1C is not detected down to the limiting sen- 
sitivity level (3cr) of ~ 7.5 x 10~ 5 Jy. Actually, the only source 
clearly detected in our L'-band image is TMR-1 AB itself. There 
is also a slight hint of nebulous emission close to TMR-1AB and 
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Filament 



TMR-1C 



Fig. 1. a) ISAAC Ks (2.16/im) image of TMR-1 and its surrounding. The total integration time in this image is 20 minutes and 
the image size is approximately 2.4' x 2.4' b) Zoom-in on the Ks-image of TMR-1 clearly showing the candidate very low-mass 
object TMR-1C at the end of the filament structure arising from the nebulosity around the binary TMR-1 AB, and our detection of 
a 'counter-object', TMR-1D. Image size is approximately 40" x 40" c) H-band image of the same field as shown in b). d) L'-band 
image of TMR-1, at the same scale and orientation as the fields shown in b) and c). 

Table 1. Infrared photometry of TMR-1 C 



Filter 


Central 


Magnitude 


Date of observation 




Wavelength 


or Flux 


(UT) 




w 


(mag or Jy) 




H 


1.65 


19.21+0.05 


Oct 01 + 02, 2000 


Kv 


2.16 


17.63±0.1 


Dec 05, 1998 


K v 


2.16 


17.53+0.05 


Oct 01, 2000 


L' 


3.78 


< 7.5 x 10" 5 


Sept 20 + Oct 02, 2000 


IRAC CHI 


3.6 


<0.12x 10~ 3 


Mar 07, 2004 


IRAC CH2 


4.5 


<0.51x 10" 3 


Mar 07, 2004 



in the direction of the filament, as well as a faint detection of dif- 
fuse emission from the north-western nebulosity structure. The 
object TMR-1D was not detected. 

We have also searched the Spitzer archive for data of TMR- 
1, and found that observations have been obtained as part of 



the Spitzer/IRAC programmes No. 37 (PI: Fazio) and No. 3584 
(PI: Padgett) devoted to IRAC mappings of the Taurus molecular 
cloud. Since the spatial resolution of IRAC at its 4 wavelengths 
channels, chl (3.6mic), ch2 (4.5mic), ch3 (5.8mic), ch4 (8.0mic) 
is on the order of l"to 2"and TMR-1C is almost 10"away from 
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Fig. 2. Spitzer/IRAC image at 3.6/mi (left) and 4.5/mi (right) showing TMR-1 and its surrounding. North is up and East is to the 
left, and the image sizes are the same as in FigureQJ). 



TMR-1 AB, TMR-1 C could have been detected if it is bright 
enough at IRAC wavelengths. Therefore, we retrieved the post- 
BCD data from the Spitzer archive, reduced with the Spitzer 
pipeline version S14.0, and investigated the images. The diffuse 
nebulosities of TMR-1 are well detected in the IRAC 3.6mic 
and 4.5mic images, which we show in Figure [2] and they are 
very similar in appearance as in our ISAAC 2.2mic image, de- 
spite at much lower spatial resolution. The IRAC channel 1 has 
a central wavelength close to the ISAAC L' filter, but due to its 
much larger pixels IRAC is more sensitive to extended emission 
than ISAAC. At IRAC channel 3 and 4 (5.8mic and 8.0mic) only 
the main source TMR-1 AB is seen. On none of the mosaiced 
IRAC images TMR-1 C is detected. This may partly be due to 
the closeness of TMR-1 C to the very bright main source TMR- 
1AB, which has a flux of ~1 lOmJy at 3.6mic, and ~145 mJy at 
4.5mic, but may also be related to an intrinsic faintness of TMR- 
1C in this wavelength range. Given the local background and the 
Poisson noise measurement at the expected position of TMR-1 C, 
we estimate upper 3<x flux detection limits of 0.15mJy (3.6yum) 
and0.51mJy (4.5/im). Concerning the 5.8yum and 8. Oyum images, 
the PSF wings of the extremely bright source TMR-1 AB prevent 
any meaningful determination of an upper flux limit. The upper 
flux limits extracted from the Spitzer data are listed together with 
the ISAAC photometry in Table Q] 



2.3. VLT/ISAAC K-band Spectroscopy 

K-band low -resolution (R=450) long-slit spectroscopy of TMR- 
1C, and of the filament was performed with ISAAC during the 
night Dec 04-05, 1998. In Figure [3] we show the slit position- 
ing on TMR-1 C and the filament. The width of the slit was 
l"and spectra were obtained by observing the object alternately 
at two positions along the slit, i.e. by nodding along the slit with 
a nodthrow size of 40". An exposure time of 120 seconds per 
individual frame was used and the resulting effective exposure 
time of the final co-added spectrum is 7680 seconds. Since the 
slit was positioned in a way to include TMR-1 C and most of 
the filament, the binary TMR-1 AB was correspondingly slightly 




Fig. 3. Position of the slit on TMR-1 C and the filament is shown. 
The sizes and positions of the apertures for the extracted spectra 
are indicated. The underlying image is the ISAAC Ks aquisition 
image on logarithmic scale. 



out of the slit and only light from its PSF wings and from the 
immediate circumbinary environment entered the slit. 

The spectral images were sky-subtracted in pairs using the 
respectively nodded companion image for subtraction, cleaned 
of pixels affected by cosmic rays and non-linear response, and 
then combined to a single spectral image of TMR-1, which is 
shown in Cuby et al. (2000). Four different apertures have been 
extracted from this 2D-spectrum at positions that are shown in 
Figure [3] The largest aperture (H2 knot 1) actually includes two 
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separate peaks of molecular hydrogen emission (cf. Figure [8]). 
Both knots appeared to be very similar on the 2D spectral image 
and a single aperture was therefore used to extract the spectrum, 
improving the signal to noise ratio. 

In order to correct for telluric absorption, as well as to at- 
tempt flux-calibration, the spectrum of a B6IV standard star 
(BS3672) was observed directly after TMR-1, with the same in- 
strumental setting and in the same manner (nodding along the 
slit) as the science spectra. After the basic reduction the stan- 
dard star spectrum was divided by a black body spectrum cor- 
responding to the spectral type of the star from what the abso- 
lute spectroscopic response was derived. The extracted TMR-1 
spectra were then divided by the spectroscopic response func- 
tion. Absolute flux calibration based on the spectroscopic obser- 
vations of the standard star had shown to be difficult, because 
it relied on the assumption that all flux from the standard star 
was inside the slit. Moreover, the seeing worsened by a factor of 
~ 1 .5 from the time of observation of TMR- 1 to the standard star 
observation. Consequently, the spectra were calibrated to match 
the flux density of TMR-1 C as determined from its photometric 
measurement (Sec. 12. II ). Wavelength calibration was made using 
the OH sky emission lines and telluric absorption lines. The esti- 
mated accuracy in wavelength calibration is on the order of one 
pixel, i.e. ~ 7A. All final four spectra are plotted in Figure[4] 

The following additional uncertainties affect the resulting 
spectra. Since the standard star spectrum has not been corrected 
for hydrogen absorption lines, these lines appear in emission in 
the TMR-1 spectra after division by the response function. In 
the case of the spectrum close to TMR-1AB and the spectrum 
H2 knot 1, however, Bry emission is also seen in the raw spec- 
trum, while no emission is found in the spectrum H2 knot 2. 
Therefore, the Bry line in the first two spectra of Figure 3 is 
overestimated. In order to roughly estimate the actual intensity 
of the Bry (2.16ji/m) emission line in these spectra relative to 
the 2.12yum emission line, the response function was corrected 
for Bry absorption through a simple interpolation over the re- 
gion /l[2.15/mi,2.175yum]. This way we estimate that the actual 
intensity ratio between the Bry and the 2.12/mi emission line is 
~ 0.7 in the spectrum close to TMR-1AB and ~ 0.14 and the 
H2 knot 1 spectrum. Absorption of the Earth's atmosphere does 
significantly affect the quality of the spectra at and below 2/im, 
and also to some extent above 2.4/im. For instance, below 2/vm 
the spectra are noisier by factors of 4-18 compared to the spec- 
tral region A > 2fim. 

Fortunately, a large number of emission lines that we iden- 
tified in the TMR-1 spectra arise from ro-vibrational transitions 
of molecular hydrogen, which occupy regions of the spectra that 
are not affected by telluric absorption nor by hydrogen absorp- 
tion lines of the standard star. Therefore the spectral signatures 
are well measurable. The relative line fluxes of identified emis- 
sion lines, normalized to the H2 1-0 S(l) flux, have been mea- 
sured and are listed in Table|2] 



3. Results and discussion 

3.1. Morphology of the nebular structure, new features and 
detections 

It is well known from previous observations that TMR-1 is asso- 
ciated with bright extended nebulosity consisting of mainly two, 
a north-western and a south-eastern structure. Our ISAAC im- 
ages show these structures in great detail (FigureQJ), c). Each of 
the two main nebulosities show a filament arm extension, with 
the filament arm of the north-western nebulosity being more dif- 



fuse than the sharp filament of the south-eastern structure. At the 
end of each filament arm we detect respectively a point-source, 
TMR-1C and TMR-1D. Another patch of diffuse emission, with 
an almost rectangular shape, is seen in the Ks-image, located 
~ 10" north of TMR-1AB. It is significantly less bright in the 
H-image, suggesting that this structure is more highly extincted. 

Considering that TMR-1 is a young stellar object during an 
early evolutionary state, the extended near-infrared emission is 
interpreted as scattered light from circumstellar dust agglomer- 
ations within a dusty proto-stellar envelope and from patches 
of molecular cloud material remaining from the proto-stellar 
collapse. Indeed, using near-infrared polarimetric observations 
Whitney et al. (1997) report a high degree (~50%) of polarized 
K-band emission from both of the main nebulosities and parts of 
the filaments, indicating that light from the central source TMR- 
1 AB is scattered off from them, and thereby producing most of 
the emission detected in the H- and K-band images. The observa- 
tions imply that the material of the NW and SE nebulosity must 
be physically close to TMR-1 AB, because otherwise they would 
not appear in scattered light. A bit surprisingly, the NW nebula 
appears much fainter at 1.65yum than at 2.16/im and also slightly 
redder than the SE nebula. We determine H-K « 2.6 mag for the 
NW nebula and H-K « 2.45 mag for the SE nebula (dominated 
by the source TMR-1AB). A large clump of dense gas that sits 
partly in front of the NW nebulosity, detected in the CS molec- 
ular line observations of Ohashi et al. (1996), is the most likely 
cause for the additional reddening. 

However, the overall picture of TMR-l's protostellar enve- 
lope+disk+outflow structure is quite complex and appears un- 
determined, despite the availability of a number of observations 
at various wavelengths. What seems confirmed is that TMR- 
1 AB currently drives a molecular outflow that had been detected 
at different molecular lines, with an outflow position angle of 
~ 170 ± 10° (Hogerheijde et al. 1998, Terebey et al. 1990). Red- 
shifted gas is found north of TMR-1AB, and the NW nebulos- 
ity appears associated with lower velocity red-shifted gas. Thus, 
the NW nebulosity and its filament has been interpreted as part 
of an outflow cavity. A cavity opening angle of ~ 30° was de- 
rived from multi-wavelength envelope modeling by Furlan et al. 
(2008). There is no clear detection of a blue-lobe of the outflow, 
which would be expected south of TMR-1AB. Instead, most 
of the blue-shifted gas seems peaked at or close to TMR-1AB, 
that is roughly at the center of the SE nebulosity. The SE nebu- 
losity, therefore, has been suggested to trace the gravitationally 
bound circumstellar or circumbinary material close to the proto- 
stars. The filament pointing towards TMR- 1C, on the other hand, 
could be part of an outflow cavity rim. 

Adding to the complexity, we clearly observe in our deep 
ISAAC images a striking apparent symmetry between the NW 
and SE nebulosities, the extending filament arms and the pres- 
ence of a point-source on either end of the filaments (cf. 
Figure QJ), c). The center of symmetry is determined by a large 
"gap" or lane of 4'.' 5 in size, i.e. almost 700 AU, in between 
the two nebulosities (see also, Terebey et al. 1990, Itoh et al. 
1999). While in some previous studies the dark lane had been 
interpreted as an extinction structure, we suggest that the gap is 
rather devoid of high density material, because neither cold dust 
emission maps from either single dish or interferometer observa- 
tions nor other high column density tracers detect the dark lane 
structure (Motte & Andre 2001, Petr-Gotzens et al. 2002, Walsh 
priv. comm.). The same conclusion was reached by Whitney et 
al. (1997) based on the relatively blue H-K colour observed to- 
wards the gap structure. TMR-1AB appears offset from the gap. 
This may be a geometrical effect, if TMR-1AB would be sur- 



Petr-Gotzens et al.: The circumstellar environment of the YSO TMR-1 and the candidate very low-mass object TMR-1C 




1.« 



2.0 2.1 2.2 2.3 

Wavelength (fim) 



2.4 



2.5 



25 



20 



2 10 



H, knot 1 




1.9 2.0 2.1 2.2 

Wavelength (//m) 



2.3 



2.4 



2.5 




l.i 



2.0 2.1 2.2 

Wavelength (/i.m) 



2.3 



2.4 



2.5 




2.0 2.1 2.2 

Wavelength (fj,m) 



2.3 



2.4 



2.5 



Fig. 4. Extracted TMR- 1 spectra. All spectra cover the same wavelength range from 1 .85/im-2.5/im. The identification and positions 
of ro-vibrational H2 lines and the Bry emission line are indicated in the spectrum of the H2 knot 1 . Note that atomic hydrogen 
emission lines (in particular the Bry-line) are caused or reinforced by division through the standard star, which shows hydrogen 
absorption (see text for further discussion). 
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Table 2. Observed line fluxes in TMR-1 spectra, normalized to the H2 1-0 S(l) flux. 



Spectral 


Measured 


Normalized 


line 


wavelength 
(A) 


integrated line flux 


H 2 1-0S(1) 


21222 


1.00+0.02 


Bry (7-4) 


21664 




close to H 2 1-0S(0) 


22231 


0.22 ±0.03 


TMR-1AB H, 1-0Q(1) 


24057 


1.29+0.29 


Ho 1-0 Q(3) 


24226 


0.86 ±0.11 


H 2 1-0 Q(4) 


24377 


0.33 ±0.06 


H 2 1-0 Q(5) 


24536 


0.69 ±0.09 


H 2 1-0 S(3) 


19573 


0.49 ±0.01 


H 2 1-0 S(2) 


20332 


0.34 ±0.02 


H 2 1-0 S(l) 


21211 


1.00 ±0.004 


Bry (7-4) 


21659 




H 2 1-0 S(0) 


22220 


0.23 ±0.01 


H 2 knot 1 H 2 2-1 S(l) 


22466 


0.12 ±0.002 


H 2 1-0Q(1) 


24044 


1.06 ±0.02 


H 2 1-0 Q(2) 


24110 


0.29 ±0.01 


H 2 1-0 Q(3) 


24214 


0.87 ±0.01 


H 2 1-0 Q(4) 


24353 


0.46 ±0.02 


Ho 1-0 Q(5) 


24524 


0.80 ±0.01 


H 2 1-0 Q(6) 


24724 


0.40 ±0.01 


H 2 1-0 S(3) 


19571 


0.71 ±0.03 


H 2 1-0 S(2) 


20328 


0.26 ±0.05 


H 2 1-0 S(l) 


21206 


1.00 ±0.04 


H 2 1-0 S(0) 


22221 


0.20 ±0.01 


H 2 knot2 H 2 2-1S(1) 


22465 


0.18 ±0.10 


Hz 1-0 Q(l) 


24039 


1.22 ±0.02 


H 2 1-0 Q(2) 


24099 


0.30 ±0.07 


H 2 1-0 Q(3) 


24208 


1.14 ±0.23 


H 2 1-0 Q(4) 


24351 


0.58 ±0.13 


H, 1-0 Q(5) 


24522 


0.85 ±0.11 



Note: No line fluxes are given for the Bry line, because the observed fluxes are overestimated (see main text for further discussion). 



rounded, for example, by a large ring of dense gas with an inner 
hole and the ring being tilted towards the observer. The nebu- 
losities would trace to some extent the optically thick surface of 
such a ring. But without more detailed information on the gas 
kinematics and observations at multi-wavelengths, no firm con- 
clusion can be made. 

In the context of our study of the nature of the point-source 
TMR-1 C, yet the most noteworthy observation revealed by our 
images is the symmetric location of TMR-1 C and TMR-1D, 
each at the end of a long filament structure. We interpret the 
morphology as such as if the filament arms had been swept up 
by a common symmetric outflow. The filaments are not an out- 
flow or jet by themselves, because they are composed of large 
amounts of continuum emission. Maybe the expansion of a pow- 
erful bipolar outflow had triggered the collapse of small low- 
mass cores, leading to the formation of TMR-1C and ID? Based 
on our ISAAC images, we measure a position angle of ~ 145° 
for a line tracing roughly the direction of the filaments and con- 
necting TMR-1 C with TMR-1D. The position angle of the cur- 
rently observed molecular outflow, detected at several different 
molecular lines has been determined to ~ 170± 10° (Hogerheijde 
et al. 1998, Terebey et al. 1990). One may speculate that the main 
outflow direction has changed counter-clock wise over the past, 
because the current outflow's axis seems different from the po- 
sition angle formed by TMR-1 C and ID. Or alternatively, one 
binary component's outflow has ceased already and we are cur- 



rently only witnessing the other component's outflow. In this 
context we note another intriguing new detection from the deep 
ISAAC image: clearly visible in FigureQJ, detected at 1.2' dis- 
tance south-southeast of TMR-1, there is a bright arc-like struc- 
ture of likely swept-up material, whose morphology resembles 
that of a bow-shock. However, no excessive line emission at 
2.12yum is detected, which should be expected if the emission 
would be caused by shock-excitation. The SpitzerflRAC images 
at 3.6mic and 4.5mic confirm the structure (Figure [2]), but nei- 
ther the IRAC images nor our deep ISAAC images show an as- 
sociated north-northwest counterpart. This is unlikely due to in- 
creased extinction, because there is no indication for this, nei- 
ther from star counts nor from dust continuum maps. Although, 
clumpy extinction on small scales cannot be excluded. If a coun- 
terpart feature exists, its emission must be faint. The arc struc- 
ture may represent a fossil, compressed post-shock region, i.e. a 
signature from an earlier passage of a shock wave. Its position 
angle, measured from TMR-1 AB to the arc head is ~ 160°. In 
Sec. |4] we will follow-up on the discussion of possible forma- 
tion mechanisms for TMR-1 C and TMR-1D, and for the fila- 
ment structures. Given these new findings we reconsider in the 
following sections the background star versus very low-mass ob- 
ject debate for TMR-1C. 
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Fig. 5. Upper spectrum: K7V star spectrum from the stel- 
lar spectral library of Pickles (1998), smoothed and extincted 
by A v =23mag, and noise added. Middle spectrum: TMR-1C. 
Lower spectrum: A T e ff = 1900K synthetic spectrum with 
log g=3.5 from the DUSTY model (Chabrier et al. 2000, Allard 
et al. 2001). Spectra are in arbitrary flux units, with the K7 V 
spectrum and the DUSTY model spectrum being offset for com- 
parison. 



3.2. TMR-1 C, a background star ... ? 

Since TMR-1 is embedded in a region of the Taurus molecular 
cloud which is heavily obscured by dust, we test the hypothesis 
that TMR-1 C is an unrelated star, located behind the molecu- 
lar cloud and having its light extincted by A,, > 20 mag. Low- 
resolution spectroscopic observations of Terebey et al. (2000) 
actually showed that the near-infrared spectrum of TMR-1 C is 
not inconsistent with an extincted background dwarf star spec- 
trum. 

Our ISAAC spectrum of TMR-1 C is shown in Figure g] 
Its signal-to-noise ratio is about 8-10. The interpretation is 
not straightforward, as the spectrum lacks any spectral features 
which could be useful for diagnostics. The emission line seen at 
2.12/mi, and those at wavelengths beyond 2.4/im are residuals 
from molecular emission from the tip of the filament (H2knot2) 
and are not originating from TMR-1 C itself. The ISAAC spec- 
trum of TMR-1 C is very similar to the spectrum obtained at 
Keck by Terebey et al. (2000), and hence does not provide any 
additional new weight to the background star origin of TMR- 1 C. 

In Figure [5] we compare our TMR-1C spectrum with an ex- 
tincted K7V star spectrum from the spectral library of Pickles 
(Pickles 1998). The spectral resolution of the library spectra was 
degraded and smoothed to fit the resolution of the ISAAC spec- 
trum, and noise was added to resemble the S/N of the spectrum 
of TMR-1 C. The amount of extinction to be applied to the dwarf 
star spectrum was constrained by the H - K color of TMR-1C. 
With H - K - 1.65 + 0.1, as implied by our ISAAC photome- 
try and from measurements reported in T98, A v was constrained 
to 20-27 mag for normal interstellar dust extinction of the stel- 
lar photosphere's of B to early-M type dwarf stars. Applying 
an extinction equivalent of Av=23mag to the K7V library star 
spectrum, leads to a reasonable conformity with the spectrum 
of TMR-1C as shown in Figure [5] But note, that the compari- 
son was mainly based on the slope of the spectra, because no 
spectral features are present in TMR-lC's spectrum. As normal 
stars earlier than M-type don't show a significant difference in 




Fig. 6. Observed spectral energy distribution of TMR-1 C (filled 
black circles, based on Table Q] and on values reported in T98; 
filled black circles with down-heading arrows indicate upper 3<x 
flux limits). The black solid line represents the spectral energy 
distribution of a NEXTGen stellar model source (Hauschildt et 
al. 1999) with T eff =4000K, log g=5.0 and extincted by A v = 
23.0mag. The shaded band indicates the range in flux distribu- 
tion expected for objects of 3-7Mj up at the age of 1 Myr and seen 
behind Ay =10-22 mag of extinction, that are overall consistent 
with the observed SED of TMR-1C for A < 2.5/L<m. 



their spectral slope over the wavelength range analysed here, the 
spectral type of TMR-1C, assuming a background dwarf, can 
only be classified as M5 or hotter, and the K7V star spectrum 
that we show in Figure [5]just gives an example of a conceivable 
fit. 

More information on the nature of TMR-1C can be obtained 
from an analysis of its full spectral energy distribution, as much 
as it is accessible. To do so, we use our photometric measure- 
ments and photometry from T98 in order to compile the spec- 
tral energy distribution of TMR-1 C. Then, we compare to model 
spectra using NextGen stellar atmosphere models (Hauschildt et 
al. 1999), and find that for wavelengths shorter than 2.5fim, the 
flux distribution of TMR-1C is, for example, well represented 
by a main-sequence dwarf star with T eff = 4000 K, seen through 
A,, = 23mag of interstellar extinction (Figure |6j. This is con- 
sistent with the result shown in Figure [5] As mentioned already 
above, combinations of lower values of A,, and T e ff, or higher 
values of A v with higher T e ff do also lead to a good fit. With 
an upper limit of A,, = 30mag a stellar photosphere with T e ff 
as high as 20000 K is still in agreement with the observed spec- 
tral energy distribution in the near-infrared regime. However, the 
obviously strong decline of TMR-lC's flux distribution in the 
mid-infrared, as witnessed by its non-detection, is inconsistent 
with any reddened normal stellar photospheres. We conclude 
that TMR-1C is not an extincted background star. 

Our conclusion is also supported from a purely statistical 
point of view, because the chances of TMR-1C being a back- 
ground object projected along the line of sight close to TMR- 
1AB are low: we calculate the probability that a background 
object coincidently falls within 10" to TMR-1AB, which is the 
separation between TMR-1C and 1AB. For doing so, we first 
determine the stellar density of K-band sources from our rel- 
atively large field ISAAC observations, assuming that all faint 
objects are background sources. Since we detect 25 stars with a 
Ks-magnitude similar or fainter than that of TMR1C, we derive 
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~ 1 x 10 3 stars/arcsec 2 , and hence a probability of 0.22 to find a 
background object as close as 10" to TMR-1AB. But TMR-1C 
is not only found somewhere at 10" distance around TMR-1 AB, 
but at the tip of a narrow filament. We note that the 2D spectrum 
clearly reveals that the molecular hydrogen emission from the 
filament exactly stops where TMR-1 C is located. In other words, 
the statistical evaluation of the hypothesis that TMR-1 C is an un- 
related background object, needs also to add the probability for 
finding the background object at a specific direction (i.e. in the 
direction of the filament) around TMR-1 AB, which we divide 
in segments of 1 degree. Then, this probability to find a random 

o 

background star as close as 10" to TMR-1 AB and within 1 of 
a certain position angle drops to 6 x 10" 4 . 



3.3. ... or a very low-mass object? 

If TMR-1 C is not a chance projected background star, what kind 
of object is it? Since the spectrum is obviously pure continuum it 
is not consistent with a bullet nor dense knot of shocked molec- 
ular gas. Moreover, at all imaging observations reported so far, 
TMR-1C appears as an unresolved point-source; the HST image 
even shows the first Airy diffraction ring. 

The study of T98 proposed that if TMR-1 C is physically as- 
sociated with TMR-1 AB it should be a young sub-stellar ob- 
ject, since it is several magnitudes less luminous than TMR- 
1AB, which is a proto-stellar system of ~ 0.5Mq (Kenyon et 
al. 1993, T98). TMR-lC's bolometric luminosity could be as 
low as 10 3 - 10~ 4 L Q as estimated by T98 and Terebey et al. 
(2000). In Figure|6]we plot, as a shaded band, the range of SEDs 
for planetary mass objects with 3-7 Mj up and of the age 1 Myr 
(based on DUSTY and COND sub-stellar evolutionary models 
of Chabrier et al. 2000, Baraffe et al. 2003) that are consistent 
with the near-infrared spectral energy distribution of TMR-1 C. 
Depending on the mass of the object, and hence on the effec- 
tive temperature, different values of extinction, ranging from 
Av=10-22mag had to be added to achieve a fit with the near- 
infrared part of the observed SED of TMR-1 C. While the gen- 
eral slope of TMR-lC's near-infrared spectral energy distribu- 
tion is consistent with extincted photospheres of 3-7 Mj llp mass 
objects, the expected value of effective temperatures for such 
low-mass objects (< 2100K), however, is in strong contradic- 
tion with our observed K-band spectrum that indicates higher 
T e ff (<; 3000K). K-band spectra of cool objects with effective 
temperatures < 2500K display prominent signatures in the K- 
band, like deep methane and/or water vapor absorption bands, 
the v = 2 - vibration-rotation bands of CO at A > 2.29/j.m, 
and the Nal atomic absorption line at 2.21/im (Burrows et al. 
2001, and references therein). In general, the depth of the H2O 
and CH4 absorption increases from M-class to L- and T-class 
objects, though at some phase during L spectral class H2O ab- 
sorption is weakened by dust formation and re-radition. But in 
any case, the molecular absorption is such pronounced, espe- 
cially for the coldest very low-mass T-dwarfs, that it is easily 
detectable even in low signal-to-noise spectra (e.g. Cuby et al. 
1999). It further seems that these strong absorption bands are 
also present in young (< lMyr), low gravity, low-mass objects. 
For instance, distinct water vapor absorption bands have been 
observed in near-infrared spectra of brown dwarfs and planetary 
mass objects in the very young p Oph and Orion Trapezium clus- 
ter (Wilking, Greene & Meyer 1999; Lucas et al. 2001; Lucas 
et al. 2006). None of the features described above are present 
in our spectrum of TMR-1 C. In Figure [5] we show for com- 
parison the synthetic spectrum of a young 5 Mj up object having 




Fig. 7. Spectral energy distribution of TMR-1 C (indicated by 
filled dots, triangles represent upper limits) compared to a model 
SED for a central source of 0. 1 3M with a circumstellar disk in- 
clined to our line of sight by i - 87 deg, an outer disk radius of 
94 AU, an age of « 4 x 10 5 yr, and located at a distance of 138 pc 
seen through an interstellar extinction of Ay=18.0mag. 



T e ff = 1900 K and log g=3.5. This model spectrum, generated by 
the Lyon group based on the Ames-DUSTY stellar evolutionary 
model (Chabrier et al. 2000, Allard et al. 2001), is clearly dif- 
ferent from the spectrum of TMR-1 C. Furthermore, as implied 
by Figure[6] any young planetary mass object should be brighter 
than our measured upper detection limits at mid-infrared wave- 
lengths. 

An alternative solution may be that TMR-1 C is physically 
associated with TMR-1 AB, but instead of being an extremely 
low-mass object, it may be a slightly higher mass object hav- 
ing its light significantly dimmed by a surrounding circumstellar 
disk seen edge-on. In this case, the emergent flux in the near-IR 
is dominated by scattered light from the disk surface. TMR-1 C, 
if occulted by a disk, would appear subluminous for its spectral 
type and could exhibit an observed near-IR spectral energy dis- 
tribution different from the object's intrinsic colours. This hy- 
pothesis for the nature of TMR-1 C was actually already men- 
tioned in an earlier paper by Hartmann et al. (1999). Direct and 
indirect evidence for the presence of disks around young substel- 
lar objects have been found by several authors (e.g. Natta et al. 
2002, Luhman et al. 2007). Such studies show that disks around 
young brown dwarfs have properties which are, in general, sim- 
ilar to those derived for circumstellar disks of low-mass stars. 
Also, low-resolution near-IR spectra of objects with edge-on 
disks often appear featureless (e.g. Brandner et al. 2000, Scholz 
et al. 2008), as it is the case for the TMR-1 C spectrum. To ex- 
plore this possibility, we used the SED model fitting toofj by 
Robitaille et al. (2007) that has been developed to analyze the 
SEDs of young stellar objects. A particular advantage of this 
tool is that it allows the inclusion of upper flux limits. In order to 
limit the solutions to models that are physically relevant for the 
case of TMR-1C, we postulate that the source's distance must 
be within 117-157pc, the range in distance considered for the 
Taurus molecular cloud (Torres et al. 2007), and that its age must 
be <1 Myr. The model that then agrees best with the observed 
SED of TMR-1C is shown in Figure [7] The model SED shows 
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a steep decline beyond 2.2/im, and a strong excess at A >10mic. 
The underlying model is that of a central source of ~ 0.1 3M 
at a distance of 138 pc, and with the following parameters: age 
of the source ~ 4 x 10 5 yr, A v = 18 mag, T e ff=3000K, inclina- 
tion angle of the circumstellar disk, i - 87°, inner disk radius, 
r, = 1.8 AU, and disk outer radius, r - 94 AU. 

The predicted size of the disk seems large, given that TMR- 
1C appears unresolved in all the HST and ISAAC images. An 
edge-on disk as large as r - 94 AU and its associated bipolar 
scattered light structures should have appeared resolved in the 
high spatial resolution HST images (FWHM- a.' 15 - 072) and 
possibly also in the best ISAAC images. The fact that TMR-1C 
appears as a point-source implies that the actual disk must be 
smaller than ~20 AU. A decreasing outer disk size will mostly 
affect TMR-lC's spectral energy distribution longward of 10jt/m, 
i.e. in a regime that is unconstrained due to the lack of observa- 
tions. In a picture where TMR-1C had been ejected in a dynam- 
ical interaction within the small proto-stellar cluster of TMR-1, 
it would be expected that any primordial disk of TMR-1 C has 
been truncated (cf. Sec. @). We note that if one allows models 
for sources at slightly larger distances (up to ~170pc) or slightly 
older objects (up to ~ 1.2 x 10 6 yr), one finds that also objects 
with ~ 0.2 - O.4M having disk sizes as small as 19 AU provide 
reasonable fits to the SED of TMR-1 C. But since the SED is 
unconstrained beyond 5/mi, it is currently not possible to distin- 
guish any further between the model solutions. SED models for 
objects with central source masses < O.1M are not available. 
However, an important point we wish to emphasize, is that all 
models that provide a good fit imply a circumstellar disk with 
high inclination angle of i > 80°. 

Observational evidence for an edge-on disk surrounding 
TMR-1 C would be an infrared excess emission at A > lOfim. 
Unfortunately, the spatial resolution of available observations at 
mid to far-infrared wavelengths are too low to disentangle a po- 
tential emission from TMR-1 C from the emission of TMR-1 AB. 
Other evidence, as discussed above, could be provided by near- 
infrared images with high sensitivity and a spatial resolution bet- 
ter than 071 that should reveal any resolved bipolar scattered 
light structure. 

3.4. On the nature of the filament 

It has been suggested that the arc-shaped filamentary structure 
that seems to 'connect' TMR-1 C with TMR-1 AB is a material 
tail formed during an encounter of proto-stellar disks surround- 
ing TMR-1 A and B respectively (Lin et al. 1998, T98). The same 
encounter may have caused the formation of TMR-1 C from frag- 
mentation of a part of the filament into a very low-mass object, or 
through ejection during the dynamical interaction. As such there 
would be a clear physical relationship between the filament and 
TMR-1 C. 

The higher spectral resolution of our ISAAC data as com- 
pared to the previous spectrum of the filament (Terebey et al. 
2000) allows us to extract new physical information for the na- 
ture of the filament. Previous studies mainly report on scattered 
light emission arising from the filament. However, in spite of a 
clear component of continuum emission, we believe that there is 
another distinct process present, that is shocks. In Figure [8] we 
plot the spatial distribution of the intensities in the H2 1-0 S(l) 
and H2 1-0Q(1) emission lines, and in the continuum, along the 
filament. Obviously there are 2 locations where the emission in 
molecular hydrogen is significantly in excess. Spectra at these 2 
locations along the filament have been extracted (H2 knot 1 and 
H2 knot 2 in Figure|4] note: the H2 knot 1 location, i.e. aperture 2, 
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Fig. 8. Spatial distribution of the intensities in the H2 emission 
lines (dotted and dashed blue lines) along with the continuum 
(red solid line). The intensities have been normalized to their 
peak emission in the ISAAC TMR-1 spectrum. The apertures of 
extracted spectra are also indicated: (1) close to TMR-1 AB, (2) 
H 2 knot 1,(3) H 2 knot 2. 



actually includes 2 peaks of H2 emission which showed indistin- 
guishable spectra and where therefore combined to improve the 
signal-to-noise). Both spectra show a wealth of emission lines, 
which are predominantly ro-vibrational lines of molecular hy- 
drogen. 

Molecular hydrogen emission mainly originates from either 
of the two physical processes: shock excitation or ultraviolet flu- 
orescence. We use the intensity ratio of the transitions H2 v=l- 
0S(1) at 2.12yum and H 2 V=2-1S(1) at 2.24fim to distinguish 
between these two cases. Thermal excitation via shocks should 
be the responsible mechanism if this ratio is > 5, and gas densi- 
ties are not too high, i.e. < 10 5 cirT 3 (Shull and Beckwith 1982). 
Since we measurea2.12/mi/2.24yumratioof ~ 8.25 and of ~ 5.6, 
for the H2 knot 1 spectrum and the H2 knot 2 spectrum respec- 
tively, and we do not observe any lines from transitions of high 
vibrational levels, we conclude that the H2 emission in the fila- 
ment is due to shock excitation. Support for the shock interpre- 
tation is also gained from H2 ortho-para ratios discussed below. 

Slightly different is the situation for the spectrum extracted 
close to TMR-1AB, which does not show any detectable emis- 
sion at 2.24/im, suggesting the absence of molecular gas at 
higher temperatures (> 2000 K). The Bry recombination line 
of atomic hydrogen at 2.166^m, on the other hand, is one of 
the most prominent lines in this spectrum. This indicates that 
a process must be present that is capable of ionizing hydro- 
gen. Winds from the young stellar object is one of the likely 
responsible mechanisms. But also disk accretion onto the central 
source may be the cause for the Bry emission. Prato et al. (2009) 
present a K-band spectrum of TMR-1 AB which shows similar 
features as our ISAAC spectrum obtained close to TMR-1 AB 
(remember that the ISAAC slit positioning is offset by almost 1 " 
from the TMR-1AB peak emission). These authors determine, 
from the Bry emission line luminosity, a mass accretion rate of 
1 .9 x lO~ 7 M /yr for TMR-1, which is quite consistent with Bry 
emission being generated by accretion. 

What causes the shocks in the filament? It has been known 
for long, that TMR-1 is associated with a large molecular out- 
flow (Terebey et al. 1990, Hogerheijde et al. 1998) observed 
e.g. at 12 CO (3-2). The direction of this outflow is in very good 
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agreement with the direction of a jet indicated in Fell narrow- 
band near-infrared images (Petr-Gotzens et al. 2002). The po- 
sition angle of the filament is different by ~ 20° from the po- 
sition angle of the outflow and jet. A plausible scenario is that 
the filament is part of the edge of a cavity cleared by a lower 
velocity outflow, as for example observed in HH211 (Gueth & 
Guilloteau 1999). The molecular hydrogen emission knots along 
the filament would then arise, because this lower velocity flow 
hits into the cavity rim, thereby creating a C-type shock (see 
Sec. 3.4.1). Or, the filament is intrinsic, pre-existing dense mate- 
rial being shocked by the outflow as it happens to be in its way. 
It is unlikely that the filament is a jet by itself. While bent jets 
have been observed (Davis et al. 1994, Bally & Reipurth 2001), 
the significant amounts of continuum emission coming from the 
filament, together with a high degree of polarization (Whitney et 
al. 1997), and detection of high column densities associated with 
the filament (Hogerheijde et al. 1998, Motte & Andre 2001), ex- 
clude a jet. 
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3.4.1 . A shock-model for the emission from the filament 

While we could establish, that the H2 emission along the fila- 
ment is caused by shocks, the evaluation of the precise details of 
which shock excitation mechanism is at work, and what is the 
underlying shock physics is more complex. A powerful method 
to analyze different possible mechanisms is to compare mea- 
sured Column Density Ratios (CDRs) to various shock excita- 
tion models in a schematic CDR diagram. This approach was 
taken in this study, and we converted our flux measurements 
of the ro-vibrational molecular hydrogen lines (Table |2]i into 
Column Density Ratios following the calculations as given in 
Eisloffel et al. (2000) and Smith (2000). Thus, we determine 
the column of H2 in the upper energy level necessary to gen- 
erate each line flux. In the CDR diagrams, we compare column 
density ratios to those of a slab of gas of constant temperature 
(2000K). By so removing the exponential Boltzmann factor, we 
are able to analyse the results in detail. 

The CDR diagrams for each position along the filament are 
shown in Figures I9I10I11I The shock model plotted on all di- 
agrams is a C-type bow shock. The sparesness of the data, in 
particular with respect to the range in temperatures, and their 
scatter prevent detailed fits. Note that a constant temperature 
model (corresponding to a linear fit) would be consistent with 
all the data. A horizontal line on the CDR diagrams would corre- 
spond to 2000 K. A possibility is a single C-type shock in which 
the emission derives from an extended zone with a quite narrow 
range of gas temperatures (Smith & Brand 1990). 

However, shock waves of all types do not maintain gas at a 
fixed temperature and we expect a range of temperatures to be 
present. As an example, we demonstrate on the diagrams that 
at all locations a gently curved shock with some range of shock 
velocities is also consistent. This issue can be resolved by ob- 
serving pure rotational lines from low levels in the mid-infrared 
since very high CDR values are expected from bow shocks (see 
lower panels of the Figures fTOIllb . 

All CDR diagrams display no evidence for ortho-para ra- 
tios deviating from that usually associated with shock waves. 
Fluorescent emission from unshocked cloud gas in star forma- 
tion regions yields enhanced columns of gas in the para upper 
energy levels (e.g. Li et al. 2008). Here, however, we find that, 
within the errors, there is no deviation from a thermalised ortho- 
para ratio of 3 (Smith et al. 1997), corresponding to linear fits 
on the diagrams (note the large errors associated with the noisy 
parts of the high Q-branch and high S -branch transitions). 
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Fig. 9. Log(CDR) diagram derived for the position close to 
TMR-1 AB. The solid line represents a C-type bow shock model. 



3.4.2. Extinction along the filament 

Estimates of the interstellar extinction can be directly derived 
from line flux ratios of the Q-branch at 2.4-2. 5^m, independent 
of the specific shock model. 

In Figure [9] we show the data for the filament position close 
to TMR-1 AB with 2 different extinctions (upper and lower panel 
of the figure). With Ak = Omag all the Q-branch datapoints lie 
above the S-branch datapoints, while a much better fit of all dat- 
apoints with the model is achieved for Ak = 1.8mag (lower 
panel). The uncertainty in the derived extinction is +0.3 mag. 
This value can be compared with the extinction determinations 
toward TMR-1 AB reported in Prato et al. (2009) and based on 
low-resolution K- and L-band spectra. They find Ay = (24.2 + 
l.l)mag using an optical depth measurement of the ice absorp- 
tion feature at 3.1/mi, and Ay = (20.8 + 1.8)mag based on the 
source's near-infrared colors. Both extinction values are higher 
than our extinction derived from the molecular hydrogen line ra- 
tios. This is not surprising as the different methods trace different 
lines of sight through the material causing the extinction. In par- 
ticular, the molecular hydrogen emission lines are supposedly 
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Fig. 10. Log(CDR) diagram derived for the H2 knot 1 position. 
The solid line represents a C-type bow shock model. The upper 
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includes the measurement at H2 2-1 S(l). 



zz 


\ H2knot2 


D v=1, S 


X. 






8 1 -° 




O v=1, Q - 







A v=2, S . 


CN 






Ss./ 




■ 


01 




. 


8 0.5 




- 


\ 








01 






^ A^ 









■ 


" 0.0 






- 


CP 






■ 


O 




. 


^~ 


. . . ■ . . . ■ . . . 1 . 


. . 1 . . . 1 . . . 1 . . . 



4000 8000 1 2000 
Upper Energy, Tj (K) 

Fig. 11. Log(CDR) diagram derived for the H2 knot 2 position, 
close to TMR-1 C. The solid line represents a C-type bow shock 
model. The upper panel and lower panel show the same data- 
points, with the exception that the lower panel has a larger range 
of energies, and includes the measurement at H2 2-1 S(l). 



produced in shocked regions away from the protostellar photo- 
sphere. Also, our ISAAC spectrum is slightly off (~ 1") from the 
TMR-1 AB protostar system, while Prato et al. (2009) centered 
their spectrum on TMR-1 AB. 

For the filament positions H2 knot 1 and H2 knot 2 (FigurefTOl 
and flTTi the modeling is somewhat facilitated by the presence of 
the H2 2-1 S(l) line, enabling a fit over a larger range in tempera- 
ture. The calculated extinction at the position of the H2 knot 1 is 
Ak = 1 .8mag, and hence not different from the extinction we de- 
termined closer to TMR-1 AB. At the position of the H2 knot 2, 
i.e. only ~ 1'.'3 away from TMR-1C, we measure an Ak of 
2.4mag [Ay = (21. 4+2. 7)mag], which is in good agreement with 
the values of extinction estimated towards TMR-1C (cf. Sec. 3.2 
& 3.3). Moreover, the extinction determined at the H2 knot 2 
is significantly higher than at the H2 knot 1 and close to TMR- 
1AB, indicating that the extinction increases from close to the 
main source along the filament. This finding seems to be in con- 
tradiction to the conclusion reached by T98 where they observe 



an increasing flux ratio with decreasing wavelength for TMR-1 C 
relative to the protostars TMR-1 AB and they conclude that there 
is decreasing extinction towards TMR-1C. However, TMR-1C 
could be located in front of the filament while the filament is 
indeed receding into the denser parts of the molecular cloud. 

If we consider the continuum in the spectrum taken close to 
TMR-1 AB being dominated by the continuum from the proto- 
stars, then the spectra continuum slopes (cf. Figure [U, m prin- 
ciple, imply a lower interstellar extinction towards TMR-1 C 
than towards TMR-1AB. Our SED analysis in Sec. 13.31 sug- 
gests A v = 18.0mag for TMR-1C, and Prato et al. (2009) derive 
Ay = (24.2 + l.l)mag for TMR-1AB. In addition, strong excess 
emission from accretion and the (unknown) circumstellar disk 
geometry are also contributing to the continuum spectral shape 
of TMR-1AB. Disentangling the different effects is beyond the 
scope of this paper. 
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4. Formation scenarios for TMR-1 C (and TMR-1 D) 

We now discuss different models that may qualitatively explain 
the formation of TMR-1C as a low-mass object physically asso- 
ciated with the proto-binary TMR-1 AB. In such a picture of the 
TMR-1 system, TMR-1 C and TMR-1D are both very low mass 
objects (either low-mass stars or brown dwarfs) and their sym- 
metric locations on opposite sides of the main source TMR-1 AB 
shall be taken into account when considering possible formation 
scenarios. 

A scenario worth debating is that of outflow triggered star 
formation (Sandell & Knee 2001). This mechanism naturally 
yields collapsed objects within the high-pressure region asso- 
ciated with a jet or outflow impact. Suppose a string of well- 
spaced pre-stellar cores existed in this region before one of them 
collapsed to form the binary system. Provided that the develop- 
ing proto-stellar jets are directed along the axis connecting the 
cores, the subsequent impact could trigger a collapse in adjacent 
cores simultaneously on either side of the expanding outflow. 
Once the low-mass core is forced into collapse the formation 
of a low-mass object proceeds 'as normal', i.e. its mass contin- 
ues to grow via disk accretion. This model has the advantage of 
generating symmetric collapsed objects on opposite sides at the 
locations of adjacent cores in the string. Moreover, the on-going 
outflow may even drive away extended disk and envelope ma- 
terial, generating swept-up dust ridges and shells. The filament 
would represent denser material along the outflow direction, pos- 
sibly consisting of swept-up matter, that scatters light having es- 
caped from the main source TMR-1 AB through a low density 
cavity excavated by the jet/outflow. At the densest knots along 
the filament the outflow impact gives rise to shock excited emis- 
sion. Low-mass pre-stellar cores may also be produced in the 
shock-compressed layers of swept-up dusty filaments (Padoan 
& Nordlund 2004). Sandell & Knee (2001), for example, detect 
several compact mm-sources in a curving arc of dust near the 
proto-stellar object SSV13 in the NGC1333 star-forming region. 

Intriguingly, the protostellar outflow HH211 also possesses 
two symmetric star-like objects at the tips of the presently visible 
outflow (see Fig. 1 of Eisloffel et al. 2003). As with the TMR-1 
outflow, the near-infrared H2 outflow of HH211 raises curiosity 
because it also contains a filament of material with strong con- 
tinuum emission (O'Connell et al. 2005). 

Another possibility to discuss is the formation of TMR-1 C 
(and TMR-1D) from disk fragmentation processes. Numerical 
simulations recently published by Stamatellos & Whitworth 
(2009) show that a number of very low-mass objects with masses 
in the range 5-200 Mj can be produced by gravitational disk frag- 
mentation occuring in an extended circumstellar disk. Most of 
the produced objects are brown dwarfs (>13Mj) and low-mass 
hydrogen-burning stars. Already ~ 10 4 yrs after the formation 
of these objects from disk fragmentation, dynamical N-body in- 
teractions lead to an ejection of most them, i.e. making them 
unbound and liberate them into the field, but a few objects typ- 
ically remain bound. Quite interestingly, the most common (i.e. 
the most likely) configuration resulting from the simulations af- 
ter a few hundred thousand years is that of a solar-like primary 
with a nearby low-mass stellar companion ('TMR-1AB'), and 
two wide (r>100AU) companions, consisting of a low-mass star 
('TMR-1C') and a brown dwarf ('TMR-1D'); a configuration 
that very closely resembles what we believe we observe for the 
TMR-1 system. In this model, TMR-1 C and TMR-1D remain 
weakly bound to TMR-1 AB. They are also expected to retain 
their circumstellar disks, although at relatively small sizes, as the 
disks may had been stripped off even further during dynamical 



interactions. Dynamical interactions are also held responsible for 
inhomogeneous orientations of the young objects' circumstellar 
disks. Non-coplanarity is frequently observed in the modelling 
(Stamatellos & Whitworth 2009). It seems plausible that a close 
passage of, e.g. TMR-1C to TMR-1AB, early in the evolutionary 
history of the system, caused the circumstellar disk and hence 
the jet of TMR- 1 A (or B) to precess, which is what our observa- 
tions of a bow structure detected at > 1' distance from TMR-1 
indicate. 

A shortcoming of the 'disk fragmentation with follow up 
ejection' scenario for the origin of TMR-1 C and D is that it does 
not explain the existence of symmetric filaments and a physical 
relation of those to TMR-1C and TMR-1D. Although filamen- 
tary patterns of the fragmenting primordial disk are observed in 
the simulations, they are usually short-lived and much smaller in 
size than what is observed for TMR-1. Therefore, we consider 
for a moment that both, primary and secondary component of 
TMR- 1 AB originally condensed out of the collapsing pre-stellar 
core, and that both were accompanied by their own extended 
protostellar disk. Presume that in a subsequent encounter of the 
systems, tidal structures consisting of disk material developed 
('the filaments'), and that one proto-stellar component was cap- 
tured by the more massive component to form the proto-binary 
TMR-1AB. Naturally, the tidal filaments form symmetric struc- 
tures. As frequently observed in different numerical simulations, 
very low-mass objects are then formed out of the densest parts 
of the filamentary tidal arms (Shen & Wadsley 2006, Lin et al. 
1998), often far from the main source. This way there would be a 
clear causal dependence for the filaments and the candidate low- 
mass objects TMR-1C and ID. Shen & Wadsley (2006) further 
show that the formed 'proto brown-dwarfs' posses circumstel- 
lar disks of sizes up to ~20AU. The objects have typical masses 
in the range 2 - 73Mj, but the simulations are stopped before 
mass grow through disk accretion could be completed, so that 
the objects' final masses may be slightly higher. 

5. Summary and Conclusions 

In this paper we have used VLT/ISAAC infrared imaging and 
spectroscopy data, complemented with near-infrared HST and 
Spitzer/IRAC photometry in order to study the nature of low- 
luminosity objects detected in the close surrounding of the 
young stellar object TMR-1 (actually a binary source: TMR- 
1 AB). We have also discussed the particular morphology of the 
dust distribution around TMR-1AB as it appears in the near- 
infrared images, and how the specific dust structures may be re- 
lated to a low-mass star formation process. Our main objective 
has been to add further constraints on the potential sub-stellar 
nature of the faint source TMR-1C, which is a wide compan- 
ion to the young proto-stellar binary TMR-1AB. The striking 
location of TMR-1C at the tip of a narrow filament emanating 
from the nebulosity surrounding TMR-1AB suggested a physi- 
cal association and implied a substellar nature for TMR-1C in 
previous studies (T98). Our most important results derived from 
the analysis in this paper are summarized as follows: 

- The near-infrared K-band spectrum of TMR-1C is feature- 
less, and fully consistent with previous lower resolution 
spectra. The spectrum indicates that TMR-lC's effective 
temperature is > 3000K and rejects the possibility of a very 
cold low-mass object. 

- The near to mid-infrared spectral energy distribution of 
TMR-1C clearly shows a decline in flux for A > 2.5/im. It 
can not be fitted with an extincted background dwarf star, 
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nor with a brown-dwarf or planetary mass object photo- 
sphere, for which the object should have been detected at 
mid-infrared wavelengths at the sensitivity level of our data. 

- The shape of the spectral energy distribution of TMR-1 C 
is consistent with a model SED of a very low-mass stellar 
source of 0.1 3M Q having a circumstellar disk seen almost 
edge-on (i = 87°) and located at the distance of the Taurus 
molecular cloud. The presence of the edge-on disk that sig- 
nificantly dims the central source plus additional interstellar 
extinction of Ay ~ 18 mag is then the cause for the very low 
apparent brightness of TMR-1C. Hence, the observations at 
near-infrared mainly show TMR-1 C in scattered light. 

- Through our very sensitive H- and Ks-band observations 
with very good spatial resolution we were able to detect 
a faint point-source that we call TMR-1D (Ks~18.7mag). 
TMR-1D is located roughly 20" north-west of TMR-1AB 
and, in analogy to TMR-1 C, is seen at the end of a dif- 
fuse filament-like dust structure which arises from the north- 
western nebulosity associated with TMR-1. There is an ap- 
parent symmetry of TMR-1 C and 'its filament' with TMR- 
1D and 'its filament'. TMR-1D is a very low luminosity ob- 
ject, if it is physically associated with TMR-1. 

- The candidate low-mass objects TMR-1 C and TMR-1D may 
have been formed in an outflow triggered collapse of a string 
of low-mass cores, or from fragmentation in dense tidal fila- 
ments that were generated during an encounter between pro- 
tostellar disks of TMR-1 A and B. In these scenarios the fil- 
aments are either dust material swept-up by the outflow, or 
tidal arms formed from interacting proto-stellar disks. 

- Our analysis of spectra extracted at several locations along 
the filament that 'points' towards TMR-1C shows the pres- 
ence of strong molecular hydrogen emission due to shock ex- 
citation. The physical nature of the shocks is well described 
by slightly curved C-type bow shock models. 

- The extinction towards the filament is likely increasing from 
the filament's origin, close to the source TMR-1 AB, along 
the filament in the direction of TMR-1C. This is implied 
by extinction values derived from the Q-branch of molec- 
ular hydrogen at 2.4 - 2.5//m and means the filament struc- 
ture is probably receding into the denser parts of the molec- 
ular cloud. TMR-1 C itself, however, may be sitting slighty 
in front of the filamant's tip, because there is no clear indica- 
tion that the extinction along the line of sight to TMR-1 C is 
higher than towards TMR-1AB. We suggest the filament is 
best interpreted in terms of pre-existing dense material being 
hit by the current outflow, which causes shock emission. An 
additional component of the filament is continuum emission 
from scattered light. 

Our results provide new arguments in favour of TMR-1 C be- 
ing a young low-mass object associated with the TMR-1 proto- 
stellar system rather than being an extincted background star. 
Moreover, we argued that the faint object TMR-1D may also 
be a component of the TMR-1 system. Still, unambiguous proof 
needs to be supplied. A confirmation of an edge-on circumstellar 
disk surrounding TMR-1 C would be the detection of an infrared 
excess at A > 10/im and polarized scattered light emission from a 
disk and/or an envelope. The expected flux from a disk at around 
20/im is at a level that is detectable with currently available mid- 
IR instrumentation on 8-10m class telescopes, and follow-up 
observations are crucially demanded. An associated bipolar scat- 
tered light structure, however, is probably much more difficult to 
observe, since the potential size of the disk (and probably also 
envelope) is small and is estimated to be < 0.1", based on the 



point-source morphology at all HST and ISAAC images. The 
nature of TMR-1 D is currently undetermined and follow-up ob- 
servations are also needed for this object to fully clarify its na- 
ture. 

There are currently only very few proto-stellar systems 
known to harbor very low-luminosity wide companions in their 
immediate surroundings. Such systems are especially important 
for probing very low-mass star formation models, and TMR-1 is 
potentially one of those prime-target proto-stellar systems. 
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